Domain wall creep in epitaxial ferroelectric Pb(Zr .2Tio.8)03 thin films 
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Ferroelectric switching and nanoscale domain dynamics were investigated using atomic force mi- 
croscopy on monocrystalline Pb(Zro.2Tio.s)03 thin films. Measurements of domain size versus writ- 
ing time reveal a two-step domain growth mechanism, in which initial nucleation is followed by 
radial domain wall motion perpendicular to the polarization direction. The electric field depen- 
dence of the domain wall velocity demonstrates that domain wall motion in ferroelectric thin films 
is a creep process, with the critical exponent \i close to 1. The dimensionality of the films suggests 
that disorder is at the origin of the observed creep behavior. 

PACS numbers: 77.80.-e, 77.80.Dj, 77.80.Fm 



Understanding the propagation of elastic objects 
driven by an external force in the presence of a pinning 
potential is a key to the physics of a wide range of sys- 
tems, either periodic, such as the vortex lattice in type II 
superconductors Q , charge density waves @ and Wi smer 
crystals [||, or involving propagating interfaces, such as 
growth phenomena [Q , fluid invasion || or magnetic do- 
main walls In particular, the response to a small 
external force is of special theoretical and practical in- 
terest. It was initially believed that thermal activation 
above the pinning barriers should lead to a linear re- 
sponse at finite temperature Q. However, it was subse- 
quently realized that a pinning potential, either periodic 
[|l| or disordered, jj], ||, [|, [To|, can lead to diverging bar- 
riers and thus to a non-linear response, nicknamed creep, 
where the velocity is of the form v oc exp(—/3R(f c /f)^). 
f3 is the inverse temperature, R a characteristic energy 
and f c a critical force. The dynamical exponent \i re- 
flects the nature of the system and of the pinning po- 
tential. Despite extensive studies of the creep process 
in periodic vortex systems [0, precise determination of 
the exponents has proven difficult, given the many scales 
present in this problem and the range of voltage needed 
to check the creep law Q. For interfaces, a quantitative 
check of the creep law has been done recently in ultra- 
thin magnetic films ||, where an exponent \i = 0.25 has 
been measured in very good agreement with the expected 
theoretical value for this system. Quantitative studies of 
this phenomena in other microscopic systems with other 
pinning potentials are clearly needed. 

In this respect, ferroelectric materials are of special in- 
terest. These systems possess two symmetrically equiva- 
lent ground states separated by an energy barrier [Jo, as 
illustrated in Fig. |l| for a tetragonal perovskite structure. 
Each state is characterized by a stable remanent polar- 
ization, reversible under an electric field. Regions of dif- 
ferent polarization are separated by elastic domain walls. 
The application of an electric field favors one polarization 
state over the other, by reducing the energy necessary to 
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FIG. 1: Schematic of a perovskite ferroelectric, characterized 
by two oppositely polarized ground states, separated by an 
energy barrier Uo. For Pb(Zro.2Tio.s)03, the corner Pb ions 
and the center Ti/Zr ions are positively charged, and the face 
O ions are negatively charged. 



create a nucleus with a polarization parallel to the field, 
and thus promotes domain wall motion. In addition to 
pure theoretical interest, understanding the basic mech- 
anism of domain wall motion in ferroelectrics has prac- 
tical implications for technological applications, such as 
high-density memories. In bulk ferroelectrics, switching 
and domain growth were inferred to occur by stochas- 
tic nucleation of new domains at the domain boundary, 
a behavior observed in BaTi03 and triglycine sulphate 
by studies using combined optical and etching techniques 
]l2| , |l3| , Q . Early analysis of such motion reported a field 
dependence of the domain wall speed, v ~ exp [— 1/E]. 
This bulk system behavior was explained by assuming 
that domain walls propagate via nucleation along 180° 
domain boundaries. 

In this Letter, we report on studies of ferroelec- 
tric domain wall motion in atomically flat single crys- 
tal Pb(Zro.2Tio.8)03 films by atomic force microscopy 
(AFM), allowing noninvasive investigation of domain dy- 
namics with nanometer resolution. In this model sys- 
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tern, we identify domain wall motion to be a disorder- 
controlled creep process. The dynamical exponent [i is 
found to be close to 1. The activation energy increases 
significantly from 0.5MV/cm to 1.3MV/cm as the film 
thickness is reduced from 810Ato 290A. 

The ferroelectric materials investigated were epitax- 
ial c-axis oriented Pb(Zro.2Tio.8)03 thin films, RF- 
magnetron sputtered onto conducting (100) Nb-doped 
SrTiOa substrates |H| [l(|. This system allows precise 
control of film thickness and crystalline quality, exhibit- 
ing atomically smooth surfaces with a polarization vec- 
tor parallel or antiparallel to the c-axis jl5], [HJ. To 
study switching dynamics in these films, we used a con- 
ductive AFM tip to artificially modify domain structure 
[|l5| |l7| , [l8f . Domains were polarized by applying a volt- 
age pulse across the ferroelectric film, between the tip and 
the substrate. Their sizes were subsequently measured by 
piezoelectric microscopy [[l7| , as a function of pulse width 
and amplitude. For every pulse width, we used 12 V 
pulses to polarize an array of 16 domains, and calculated 
an average domain size based on their vertical and hori- 
zontal radii. The rms error was ~10%. All domains stud- 
ied were written in a uniformly polarized area. Fig. ||a 
shows the domain radius as a function of pulse width, and 
three piezoelectric images of ferroelectric domain arrays 
written with 50/is, 1ms, and 100ms voltage pulses. As 
can be seen, varying the writing time (the pulse width) 
markedly changes the size of the AFM- written domains. 
We observe that domain radius increases logarithmically 
with increasing writing time for times longer than ~20/is. 
Below 20/iS, and down to 100ns, the shortest times in- 
vestigated, domain radius is found to be constant and 
approximately equal to 20nm, as shown by the shaded 
area in Fig. ||a All of our data suggest that this 

minimum domain size is related to the typical tip size 
used for the experiments, whose nominal radius of cur- 
vature is ~20-50nm. In previous work, we have observed 
that domain size also depends linearly on writing voltage, 
above a threshold related to the coercive field (HJ. A 
detailed analysis of the data reveals only well-defined ho- 
mogeneous domains with regular spacing, as can be seen 
in the piezoelectric image of a regular 90-domain array, 
written with 1ms pulses in Fig. ||b. We note that the 
topographic image of the same area is featureless, with 
a rms roughness of ~0.2nm. Within our ~5nm resolu- 
tion we do not detect any randomly nucleated domains. 
The data thus suggest a two-step domain switching pro- 
cess in which nucleation, originating directly under the 
AFM tip, is followed by radial motion of the domain wall 
outwards, perpendicular to the direction of polarization. 

To analyze this lateral domain wall motion, we note 
that the force exerted on the wall is given by the electric 
field E. To obtain the electric field distribution, we model 
the tip as a sphere, with a radius a. The potential on 
the ferroelectric surface at a distance r from the tip is 
then <I> ~ anc l the local field across the ferroelectric 
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FIG. 2: a) Domain size increases logarithmically with pulse 
widths longer than ^20/is, and saturates for shorter times 
as indicated by the shaded area. In all cases, the domains 
are homogeneous and no random nucleation is detected, as 
shown by three piezoelectric images of domain arrays on a 
370A thick film written with 50/xs, 1ms, and 100ms voltage 
pulses. The scale is the same in all images. In b), a large 
array shows that domains are only centered where the AFM 
tip was positioned during writing. 
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, V being the applied tip bias and d the film 



thickness. This equation will allow us to relate the change 
in domain size to the local electric field near the domain 
boundary. By writing arrays with different pulse widths, 
and subsequently calculating the average domain size for 
a given time, we can extract the speed of the domain 
wall as v = r ^ 2 Zt[ tl ^ an d the corresponding electric field 
E(r) where r = (r(ti) + r(t 2 ))/2. Fig. || shows the wall 
speed as a function of the inverse field for three different 
film thicknesses. The data fits well to a creep formula 



v ~ exp - 



R 



E 



(1) 



with [i = l. The exact dynamical exponent \i is found 
to be 1.12, 1.01 and 1.21 for the 290A, 370A, and 810A 
thick films respectively, with an estimated 10% rms error 
on the field p7[ . We find the effective "activation en- 
ergy" (R/(k B T))»E to be 1.321MV/cm, 1.305MV/cm, 
and 0.506MV/cm for 290A, 370A, and 810A thick films, 
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FIG. 3: Domain wall speed as a function of the inverse ap- 
plied electric field for 290A, 370A, and 810A thick samples. 
The data fit well to v ~ exp [—k~r (^^l w ^ n M = 1? char- 
acteristic of a creep process. 



typically 2 orders of magnitude larger than the applied 
fields during the polarization process. We note that dur- 
ing AFM writing, the effective field across the ferroelec- 
tric is approximately 10 times smaller than the field E 
[p8[ . This has no effect on the exponent \i governing the 
exponential velocity dependence, but decreases the effec- 
tive "activation energy" by a factor of 10, which has been 
taken into account in the values given above. 

Let us now consider the possible microscopic origins 
of the observed creep behavior. Creep phenomena are 
a consequence of competition between the elastic energy 
of a propagating interface, tending to keep it flat, and a 
pinning potential, preventing it from simply sliding when 
submitted to an external force. The dynamical exponent 
in the creep scenario depends on both the dimensionality 
of the system, and the nature of the pinning potential. 
Although creep processes are generally associated with 
the glassy behavior of disordered systems, they can also 
be observed in a periodic potential if the dimensional- 
ity of the object is larger than or equal to 2. For thick 
films, as for bulk ferroelectrics, the domain wall is a two 
dimensional object. In this case, from free-energy con- 
siderations, and neglecting the anisotropic dipole field 
present in a ferroelectric, one would expect the exponent 
to be ji = 1 || . In PbTi03, it has been shown theoreti- 
cally that domain wall energy depends upon whether the 
wall is centered on a Pb or Ti plane [^0|, |2l[, giving rise 
to an intrinsic periodic pinning potential. One possible 
explanation could thus be that the observed creep is due 
to the motion of the two-dimensional wall in this peri- 
odic potential. Note that this scenario is a generalization 
of the nucleation model developed for bulk ferroelectrics 



[ |l4[ . In order to test this hypothesis we calculated the 
size of the critical nucleus, using the formula derived by 
Miller and Weinreich ]T^| . To estimate /*, the critical 
length along the c-axis, we used the standard remanent 
polarization, lattice parameters and dielectric constant 
values for PZT , the 169mJ/m 2 domain wall energy de- 
rived for PbTiOa pl| , and the corrected values for the 
electric field across the ferroelectric film. We find I* to 
vary, depending on the field range, between 200 A and 
500A, 600A and 1100A, and 900A and 1700A for the 
290A, 370A and 810A. The critical nucleus would thus 
need to be larger than the thickness of the system. Fur- 
thermore, the effective "activation energy" calculated for 
the nucleation model is two orders of magnitude greater 
than the 0.5-1.3MV/cm determined experimentally. Cal- 
culations directly starting from the periodic potentials 
given in |2j]] lead to similar conclusions. These results 
strongly suggest that the films are in a two dimensional 
limit, and that the nucleation model or equivalently the 
motion through a periodic potential does not adequately 
explain the experimental data. 

The creep behavior in the film thus has to be due to 
disorder in the system, and thereby to the glassy char- 
acteristics of a randomly pinned domain wall, with the 
creep exponent dependent on the nature of the disor- 
der. Defects locally modifying the ferroelectric double 
well depth Uo and giving rise to a spatially varying pin- 
ning potential would lead to a "random bond" scenario 
similar to the one for the ferromagnetic domain walls |5| . 
The exponent \i would be \i = d ~ 2 ^ 2 ^ where £ is a char- 
acteristic wandering exponent and d the dimensionality 
of the wall. For one dimensional domain walls \i = 1/4 
whereas for two dimensional ones \i ~ 0.5 — 0.6, hardly 
compatible with the data. If, however, the defects induce 
a local field, asymmetrizing the double well, or if there 
are spatial inhomogeneities in the electric field, the sit- 
uation is different. For such a "random field" scenario 
C = ht H and thus leads to /i = 1 for 1 < d < 4. 
This scenario is therefore compatible with the observed 
data. However, further study determining the wandering 
exponent would be needed to ascertain this precise point 

Finally, we note the applications of these results to 
technological developments. The fact that the domain 
walls exhibit creep motion with a relatively large expo- 
nent fx = 1 implies a strong stability of ferroelectric do- 
mains in thin films, since the induced speed of the domain 
wall becomes exponentially small as the driving force 
goes to zero. All domains studied in this work were stable 
under ambient conditions for the entire 7 day duration of 
the experiment. Furthermore, sub-/im wide line-shaped 
domains were stable up to one month j2^|. Previously, 
large, regular arrays of ferroelectric domains with den- 
sities of ^6Gbit/cm 2 have been reversibly written [ [l9[ , 
and densities up to 150Gbit/cm 2 have been extrapolated 



from sizes of individual domains [24 . This work iden 
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FIG. 4: Piezoelectrically imaged high density ferroelectric do- 
main array at 28Gbit/cm 2 , written with 500ns pulses. 



tifies the key parameters controlling domain size: the 
strength, duration, and confinement of the applied elec- 
tric field. These can be exploited to increase the infor- 
mation storage density in ferroelectric arrays. By using 
short voltage pulses, for which the resulting domain size 
is independent of the writing time, and domain-domain 
separation as small as lOnm, regular arrays with densi- 
ties of the order ~30Gbit/cm 2 could be written. Fig. || 
shows such an array written on a 370A thick sample with 
a density of 28Gbit/cm 2 . Furthermore, it has been sug- 
gested [^5| that the minimum stable domain size is given 
by the film thickness. Our data, however, show that the 
minimum stable domain size, ~40nm, is independent of 
the thickness used in this study; rather, it is given by the 
area over which the electric field is applied during the 
polarization process. Therefore, by confining the field to 
a smaller area, domain wall creep will be limited, leading 
to smaller domains and hence even higher information 
densities. 

In conclusion, our studies demonstrate that lateral do- 
main wall motion in ferroelectric thin films is a creep pro- 
cess, governed by a characteristic dynamical exponent \i 
close to 1. The dimensionality of our films suggests that 
disorder is at the origin of the observed creep behav- 
ior, which inherently explains the measured exp(— 1/E) 
dependence of the domain wall speed. The activation 
energy is found to be around lMV/cm, decreasing with 
increasing film thickness. Finally, these results suggest a 
high degree of stability for ferroelectric domains in low 
electric fields and identify the key parameters controlling 
domain size. 
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